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SUMMARY 



The von Xarman-Killikan theory of laminar boundary 
layers presented in N.A.C.A. Technical Report No. 504 is 
applied to the laminar "boundary layer about an elliptic 
cylinder on which boundary- layer and pressure-distribution 
measurements were made at the National Bureau of Standards, 
An outline of the procedure of the von Ka'rma^n-Millikan 
method is given. 

Good agreement is obtained between the calculated and 
experimental results, indicating that the method may be 
applied generally to the laminar boundary layer about any 
body provided that an experimentally determined pressure 
distribution is available. It appears that for all 
Reynolds Numbers above 24,000 the separation point for the 
elliptic cylinder should occur at a constant distance be- 
hind the point of minimum pressure, provided that the 
boundary layer does not become turbulent. 



INTRODUCTION 



From consideration of the relative order of magnitude 
of the various terms in the fundamental differential equa- 
tions of flow, Prandtl has shown that the effects of vis- 
cosity at the high Reynolds Numbers common to aeronautical 
problems are of appreciable importance only in the thin 
fluid layer next to the surface of the body, thus consid- 
erably simplifying the fundamental equations. The actual 
solution of the boundary- layer equations has been a matter 
of considerable difficulty. It was shown, however, that 
separation of the flow from the surface is to be expected 
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when the boundary layer encounters rising pressures in the 
direction of flow. 

The simplest case, that of the boundary layer along a 
flat plate with zero pressure gradient, has been success- 
fully treated by Elasius. For this case no separation was 
to be expected. The treatment of the flow about "bodies 
when pressure gradients exist has thus far not been com- 
pletely successful. .Numerous methods have been devised 
but either the prediction of the location of the separa- 
tion point has been of doubtful accuracy or the methods 
have been applicable only to limited types of pressure 
dist ribut ions. 

Experiments on an elliptic cylinder have been con- 
ducted at the national Bureau of Standards to study the 
separation phenomena and, in particular, to check the the- 
ory of Pohlhausen. (See reference 1.) The results, com- 
pared with Fohlhausen ' s approximate solution of the von 
Karnian integral equation, have demonstrated that the Fohl- 
hausen method could only be depended upon to give reliable 
results when the pressure in the outside potential flow 
was decreasing; with increasing pressure separation may 
actually occur when Fohlhausen' s method, fails to indicate 
any flow separation. 

Von Karma' n and i.'illikan have recently devised a theo- 
ry of laminar boundary layers involving separation (refer- 
ence 2) that appears to have certain advantages over pre- 
vious theories. This solution of the laminar boundary- 
layer problem can be applied more generally than former 
methods; the equation that is used to determine the sepa- 
ration point is a close approximation near the surface 
where the separation condition is applied, thus permitting 
a more accurate determination of the separation point. 

The purpose of -the present study is to investigate 
the extent to which the theory presented in reference 2 
may be expected to sive reliable results. This investiga- 
tion was accomplished by applying the theory to the bound- 
ary layer about the elliptic cylinder on which the Bureau 
of Standards tests were made. A comparison is made between 
the calculated and experimental results, especially with 
reference to the location of the separation point. 
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EXPERI MENTAL DATA 



The tests (reference l) consisted of pressure-distri- 
bution and boundary- layer measurements on an elliptic cyl- 
inder. The pressure distribution was measured by a manom- 
eter connected to orifices in the surface of the cylinder. 
A hot-wire anemometer was used to make the boundary- layer 
surveys. The elliptic cylinder on which the tests were 
made had a major axis of 11.78 inches and a minor axis of 
3.98 inches. It was placed in the air stream with the ma- 
jor axis parallel to the general flow. The tests were 
made at a Reynolds Number of approximately 24,000, based 
on the length of . the minor axis. 

OUTLINE OF THE YON KARMAIT- HI L L I EAN PROCEDURE 



The solution of the boundary-layer problem presented 
in reference 2 is divided into two parts, an outer and an 
inner solution. The outer solution is obtained by trans- 
forming the boundary- layer equation into a form analogous 
to that for the conduction of heat , the solution of which 
is well known. This part of the solution is most accurate 
in the region where the boundary- layer velocities are near- 
ly equal to those in the outside stream, i.e., the outer 
part of the boundary layer. For reasons peculiar to the 
analysis, it is convenient to join the two solutions at an 
inflection point of the boundary- layer velocity profiles. 
Thus the range of this solution is taken as the region 
from the outer part of the boundary layer to the inflec- 
tion point, or to the wall if no inflection point exists. 

The inner solution is obtained by transforming the 
boundary- layer equation by certain approximations into an 
ordinary differential equation, which is most accurate 
near the wall. The inner solution is used only when the 
boundary- lay er velocity profiles show an inflection point. 
Its range is then over the region from the wall to the in- 
flection point. 

The steps involved in computing the characteristics 
of the laminar boundary layer by this method follow. 

Definition of d line n s ignless variables.- In the fol- 
lowing computations the velocity at an infinite distance 
from the cylinder U Q is considered to be the unit of ve- 
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locity. The minor axis I of the ellipse is used as the 
unit of length. These units of velocity and length lead 
to the formation of the following dimensionless quantities 

R = , Reynolds IJumbe.r (u, kinematic viscosity). 

Cp 

cp* - > nondimensional velocity potential. 

J 0 1 

* " V /IT 

^ - — - , nondimens ional stream function. 

u* = . nondimensional velocity in the boundary 
U o layer. 



U* = ~— , nondimensional velocity just outside the 
U o "boundary layer. 

-.2 s 

z *_ — - --- , nondimensional "energy defect" in 

^o ^o the boundary layer. 

In the following computations, the asterisks are 
dropped for simplicity, and nondimensional quantities thus 
formed are. always to "be understood. 

Calculation of outer _so.lut i o k . - In order to calculate 
the o'Her solution, it is necessary to know the variation 
of the square of the outside velocity U 5 with the veloci- 
ty potential cp along the surface. Reference 1 gives ex- 
perimentally determined values of U as a function of the 
nondimensional distance s . along the surface. This rela- 
tion is plotted in figure 1. The velocity potential 

s 

<P = / Uds is determined "by a graphical integration of 
o 

this curve. Both U and cp are t.hen known as functions 
of the distance along the surface. The U 2 ,cp curve, 
shown in figure 2, is determined by choosing corresponding 
values of U and cp at several points along the surface. 

The next step in the calculation is the expression of 
the U , cp curve in terms of power series. The calcula- 
tions are greatly simplified if the U 3 ,Cp curve is arbi- 
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trarily divided into two regions and a separate power se- 
ries is used to approximate the curve in each region; that 

n 

is, for 9 1a), TJ 3 = 2 b^*- 

i=0 

m 



9 ^ qp x TJ 3 = S Pi<p 



where 9i is the arbitrarily chosen value of 9 separating 
the two regions. It was found possible in this investiga- 
tion to obtain a satisfactory approximation to the experi- 
mental U 2 , (p curve by using t erms involving only up to 
the second power of <p . It appears that the U 2 , cp curve 
should be fitted with greater care in the neighborhood of 
the separation point than elsewhere. The value of 9i 
chosen was 0*4. The following approximate expressions for 
the U 2 ,9 curve were then obtained: 

for 9 < 0.4 U 2 = 7.3809 - 8.8559 2 

for 9 0.4 U 2 = 1.393 + 0.4149 - 0.1489 2 

The approximate U ,9 curve, given by the foregoing ex- 
pressions, is also shown in figure 2. 

It is now possible to compute the energy defect z^, 
at any point in the boundary layer as a function of 9 and 
. (The subscript co denotes the outer solution.) 



(9/<!0 = {Mo* + 7i<Pi6i* - 7 3 9! 2 g 2 *} 

+ 9 |M*i - h x *) + Pihi* + 2 72 qp 1 (g 1 * + g 3 *)J 

+ (p* {b 3 (h 3 - h 2 *) + B 3 h 2 *} (1) 

for 9 - 9i 

7i = *>i - Pi y 3 = \ - P 2 

For 9 ^ 9i put p lf p a , h 0 * , hj* , h 2 * , gl *, g 2 * 
equal to zero. The factors h Q , h lf h 2 , g x , and 
g 2 are functions of y|p only, and thus do not depend on 
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the particular case. The same functions as the foregoing 

d)c >|c sfs 

but formed for the variable are h 0 , h x , h 2 

a/ cp - q>i 

g x *, and g 2 * . Equation (l) corresponds to equation (19) 

in reference 2. The terms arising from the use of powers 
of Cp greater than 2 have been deleted. 

Substituting the values of the known coefficients in 
the expression for (equation (l)) f 

• Z(0 - jl.3930 h 0 * + 2.7864 6l * + 1.3931 g 3 *j-- 

+ cp J 7. .3800 hi - 6.9650 h x * - 6.9656' gj* - 5.9556 g a *| 



+ 9 2 |-3.8 550 h 3 + 8.7070 h 2 ~ 

The functions g 0 , r, L , g 3 . h o • h i» b 2» are computed 
from the following relations: 



} (2) 



■■here F(&) is the probability integral, 

-4= f Q * d$, a tabulated function. 

/ TT o 



,(£) = ~ ->= e + e (l - Til)) 

V TT V 7 



&2 

h 0 (t) = g 0 (t) 

hx(0 = e 0 U> + 6^(6) 



These f un ctions are 
in figures 3 and 4. 



plott ed 



throughout 



their useful range 
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In the subsequent calculations, the power-series ex- 
pansions of these functions are required. 



£ 0 (5) = 0.5000 - 0.5642£ + 0.188l| - 0.0564^ 



a fc 5 

- 0.056 

gl (£) = -0. 5642£ + e 2 - 0.5642| 3 + 0.0940E 5 

g 3 (l) = 0.188lE - 0.5642E 3 + 0. 6667I 4 - 0.282l£° f (3) 

h x (£) = 0.5000 - 1.12845 + I 2 - 0.376lE 3 + 0.037?| S 
h 2 U) = 0. 5000-1. 5045E+2E 2 -1. 5045E 3 +0. 6667^ 4 -0.1505E 5 _ 

It is now possible to compute the z^.ty curves. 
This computation is performed by choosing the value of cp 
corresponding to the desired distance along the surface 
and calculating z^ for several values of \jf by equation 
(2). The outer limit of the boundary layer is chosen as 
the point at which relative energy deficiency --- = 0.01 

where z 0 = IL-. The velocity is then 0.995 the outside ve- 
locity. " 2 

Prom the z^ curves and the relations 



u = u/l - - z - 
"o 



y J R = 2 / 



V z 0 



the velocity profiles in the boundary layer are obtained. 
Because the integrand in the expression for y Jli becomes 
infinite as the lower limit of integration is approached, 
it. is necessary to evaluate the integral by a combination 
of graphic and analytic methods (reference 2, p. 10). For 
values of z near z Q , that is, for small values of , 
z is expanded in a power series in "ty , using equations 
(2) and (3), and integrated analytically from \jr = 0 to 
some small value of ty, say \Jr = 0.05. Prom that point, 
the integration may be performed graphically by means of a 
planimeter or Coradi integrap-h. 

In cases where the u,y (or z,\J0 curves have no 
inflection point, this procedure completes the solution 
for the particular value of cp chosen.. When the u,y 
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curves have an inflection point, it is necessary to deter- 
mine certain quantities from the outer solution at the in- 
flection point to aid in fixing the boundary conditions 
for the inner solution. 

The power series in ^ for is again used to de- 

termine the position of the inflection point. Since the 
inflection point, in general, occurs at small values of 
\|r, powers of higher than ^ 3 may he neglected. The 
position of the inflection point is found by equating 

2.-J. = o, and solving for . This value of V is denoted 

where the subscript j designates a joining-point 



by t 



CO J 



characteristic. The values of 



'S.Z-N and Zj are then 



CO j 



found by substituting j in the power-series expansions 

for these quantities. This procedure completes the outer 
solution when the boundary- lay er velocity profiles have an 
inflect ion point . 

Iliner _solut ion . - The inner solution is given (refer- 
ence 2, equation (29)) in the form 



e . J_ /HZ (Lif 4 



B - 




+ 2 



+ sin 



_ l 



- s m 



i + sr 




(4) 



OZ _ 

where « 0 « = ^9. 



2 o - 



= Til 
2 



The subscript i denotes the inner solution; 

.2 



stant of integration proportional to ^g— 



y=o 



B is a con- 
and z n , 



' , and cp are derived directly from the outside pressure 



distribution, which has been experimentally determined. 

The valu< 

equal to 



The value of £. is taken from the outer solution and is 



z o " z coj' 
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It is to oe noted that ~-r-5- = 0 when t- = An- 
other boundary condition is that (Jy-) = (||-~) • 
This condition gives the equation 

/ail 

z 



_. 19. 



* a-^ &r /it?- c, 



for determining B, the left-hand side of which is known 
from the outer solution. The value of B thus found is 
substituted in (4) and i ± is then computed for several 
values of the ratio -A between 0 and 1. 

It will be found, in general, that is not equal 

to In order to make the final solution continuous, 

the values of are shifted by an amount € = t i j - l^j 

that is-, for the outer solution I = l w + £.. 

With this information, the £ , \{r curve can be drawn 
continuously for both solutions. As before, the u,y 
profile is determined from this curve by the following re- 
lations 

V z 0 



?or the inner solution, this integral may be expressed 
analytically. . In this case 



V 2 V -z 0 'cp v z y 



5iN 

o 



1/4 



- 1 



s m 



1 

/l+3 : 



1 - 

_ i 



- sin 



r 3 



Fast the joining point, the integration is most easily 
carried out graphically. 

2 

Since B is pro-po rt ional to f^-) , the condition 
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for separation is that B shall equal zero. The position 
along the surface for which this relation holds is deter- 
mined "by trial and error. 

EE SUITS AND DISCUSSION 



A comparison of the calculated velocity profiles and 
experimental results at several points along the surface 
of the elliptic cylinder is given in figures 5, 6, and 7. 
The calculated Pohlhausen curves, which have "been discussed 
in reference 1, are shown for convenient comparison. The 
agreement "becomes "better forward of the separation point. 
Close to the separation point the calculations show, in 
general, too low velocities at corresponding distances 
from the surface, hut the shape of the calculated separation 
point profile is in good agreement with the observed veloc- 
ity distribution at the experimentally determined separa- 
tion point. 

The positions of the calculated and experimental sepa- 
ration points are indicated on the pressure-distribution 
diagrams shown in figure 8, The calculated separation 
point was found to occur at a distance along the surface 
s = 1.92, Actual separation was observed at s = 1.99. In 
view of the many approximations and arbitrary procedures 
involved in the analysis, this agreement is considered very 
good. 

It is to bo noted that the curves given in figures 5, 
6, and 7 are independent of the Reynolds Number, except 
insofar as Reynolds Number affects the pressure distribu- 
tion. It is, of course, assumed that the boundary layer 
does not become turbulent. An example of the effect of a 
change in the pressure distribution on the position of the 
separation point is given in figure 8, where the calculat- 
ed separation point is shown for the perfect-fluid pressure 
distribution about the cylinder. In this case the position 
of the separation point is at s = 2.38, considerably aft 
the position found with the experimental pressure distri- 
bution. Figure 8 shows that, .for both pressure distribu- 
tions, the- separation point occurs at approximately the 
same distance, s = 0.7, aft the point of minimum pres- 
sure. Since it is to be expected that the pressure distri- 
butions for higher Reynolds Numbers will be between the ex- 
perimental curve shown and the perfect f luidrpre s sur e dis- 
tribution, this result indicates that for all Reynolds Num- 
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"bers above 24,000 the separation point for the elliptic 
cylinder will occur at a distance s = 0.7, aft the point 
of minimum pressure provided, of course, that the boundary 
layer does not become turbulent. 

A general view of the flow about the elliptic cylin- 
der in the N.A.C.A. smoke tunnel is shown in figure 9. 
The effect of separation on the entire flow configuration 
should be noted. 

CONCLUSIONS 



The computed and experimental characteristics of the 
laminar boundary layer about the elliptic cylinder are in 
good agreement. This agreement indicates that the method 
may be generally applied to the laminar boundary layer 
about any type of body provided that an experimentally de- 
termined pressure distribution is available. It appears 
that for all Reynolds Numbers above 24,000 the separation 
point for the elliptic cylinder should occur at a constant 
distance behind the point of minimum pressure, provided 
that the boundary layer does not become turbulent. 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , October 10, 1935. 



EEFEEENCES 



1. Schubauer, G. B.: Air Flow in a Separating Laminar 

Boundary Layer. T.E. No. 527, N.A.C.A., 1935. 

2. von Karman, Th. , and Millikan, C. B.: On the Theory 

of Laminar Boundary Layers Involving Separation. 
T.E. No. 504, N.A.C.A., 1934. 



N.A.C.A. Technical ITote No. 544 



Fig. 1 










I 

] 

j_ 






















--— 








I 

i 


















































































































































. 


— _. 




















/ 
/ 

h 4- 











1 


— — 





, 














/ 

._/. . 

j 

i 









_ 





















r~ H 










- -I 
















1 

\ 







- 




















































— , — 







































r 

— 































\ 




























.. v . 

V 

\ 






































. J 






















I 










i 






— I 







i-H 



O CO vjt CM 

• • • • . 

i-r 

fl 'a^xocxsa eptsq.no praoTSustnrpuofl 



ID 

to 



to 



CO 

cvJ 



«3 



D 
O 

a 



E3 
o 



O 

3 

m 

•rf 



of 
a 
o 

•H 

CO 

a 
a> 

S 

•H 

-d 
C 
o 



no 



O 



N.A.C.A. Teclinical Note Ho. 544 



Fig. 




N.A.C.A. Technical Bote No. 544 



Fig. 4 




N.A.C.A. Technical Note No. 544 



Fig. 5 
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Fig. 6 
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Figure 7.- Boundary-layer velocity distribution near separation 
point. 
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Fig. 9 




